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A  one-dimensional  mathematical  model  of  the  positive  electrode  of  a  sodium-iron  chloride  battery  for 
an  isothermal,  constant-current  discharge-charge  cycle  is  presented.  Macroscopic  theory  of  porous  elec¬ 
trodes  and  concentrated  solution  theory  are  used  in  the  model  to  describe  the  transport  processes.  The 
change  in  the  solubility  of  FeCh  with  position  and  time  within  the  cell  is  included  in  the  model  by  defin¬ 
ing  an  equilibrium  constant  that  is  a  function  of  the  NaCkNaAlCU  molar  ratio.  The  concentrated  solution 
theory  for  a  three-ion  system  with  common  cation  is  extended  to  account  for  a  diffusive  flux  of  a  sparingly 
soluble  ferrous  complex.  It  is  seen  that  this  flux  is  important,  especially  at  moderate  depths  of  discharge. 
The  effect  of  the  assumed  solubility  constant  Kspyea  on  the  battery  performance  is  characterized.  When 
JCSp,Feci  is  higher  than  1 06,  its  variation  does  not  change  the  short-time  behavior  of  the  system  appreciably. 
Simulations  suggest  that  the  iron  accumulates  near  the  sodium  tetrachloroaluminate  reservoir  during 
discharge.  When  charging,  the  net  movement  is  reversed.  As  a  result  of  continuous  cycling,  simulations 
predict  that  iron  is  depleted  at  this  boundary.  For  instance,  at  the  end  of  the  fifth  cycle,  the  iron  amount 
decreases  by  -1%  near  the  reservoir. 
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1.  Introduction 

A  zebra  battery  is  a  high-temperature  secondary  battery  sys¬ 
tem,  with  significant  promise  for  high-energy  density  applications 
requiring  long  cycle  life  [1-12].  These  battery  systems  also  have 
zero  self-discharge  and  are  unaffected  by  the  ambient  tempera¬ 
ture  [7,11].  The  zebra  battery  contains  a  liquid  sodium  electrode 
and  a  |3"-alumina  solid  electrolyte,  like  the  sodium-sulfur  battery 
[1-12].  It  also  contains  a  second,  molten  salt  electrolyte,  sodium 
tetrachloroaluminate  (NaAlCl4),  and  a  porous  metal/metal  chlo¬ 
ride  electrode  [1-12].  The  |3"-alumina  solid  electrolyte  only  allows 
Na+  ions  to  pass  and  it  has  essentially  zero  electronic  conductivity 
[10-12].  The  liquid  electrolyte  connects  the  ceramic  electrolyte  to 
the  metal-chloride  electrode  for  the  rapid  transport  of  Na+  [1-12]. 
The  battery  operates  in  the  range  of  270-350  °C  since  high  temper¬ 
atures  are  needed  to  keep  the  sodium  tetrachloroaluminate  molten 
[1,4,11,12].  In  addition,  the  resistance  of  the  solid  electrolyte  is  low 
in  this  temperature  range  [4,10-12]. 

Metal  chloride  cells  are  assembled  in  the  discharged  state  by 
mixing  the  metal  powder  with  NaCl  in  the  positive  electrode  and 
adding  the  salt  electrolyte  as  a  dry  powder  to  the  mixture.  After 
heating  the  cell  to  the  operating  temperature,  it  is  then  charged 
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to  generate  the  liquid  sodium  and  metal  chloride  [5-9,12].  Iron 
chloride  and  nickel  chloride  are  the  most  common  electrodes  used 
in  these  cells  [1,3,4,12].  Here  we  treat  the  iron  chloride  electrode. 

During  discharge,  sodium  ions  are  conducted  through  the 
ceramic  electrolyte  from  the  negative  electrode  and  then 
transferred  to  the  positive  electrode  through  sodium  tetrachloroa¬ 
luminate.  Sodium  reacts  with  iron  chloride  on  the  electrode  to 
produce  sodium  chloride  and  iron.  The  battery  is  fully  discharged 
when  there  is  no  iron  chloride  left  in  the  cell  [1-14].  The  reverse  of 
this  process  occurs  during  charging.  The  overall  cell  reaction  is  (1 ): 


2Na  +  FeCl2 


DISCHARGING  _ 

< — >  2NaCl  +  Fe 

CHARGING 


(1) 


with  a  thermodynamic  cell  potential  of  2.35  V  at  250  °C 
[1,2,4,7,12,15,16]. 

Sodium  tetrachloroaluminate  is  a  mixture  of  two  binary  molten 
salts,  NaCl  and  A1C13,  and  the  apparent  concentration  ratio  of  NaCl 
to  A1C13  determines  the  solubility  of  FeCl2  in  the  electrolyte  [17]. 
NaCl-rich  melts  are  typically  used  [14,15,17],  in  part  because  it 
is  desirable  to  maintain  low  FeCl2  solubility  to  minimize  redis¬ 
tribution  of  active  material.  Nevertheless,  the  iron  chloride  is 
sparingly  soluble,  and  with  increased  cycling,  it  does  redistribute. 
The  migration  of  the  metal  in  the  cell  results  in  a  loss  in  the  per¬ 
formance  of  the  battery,  and  may  be  a  crucial  failure  mechanism 
[9,12,15,18,19]. 

In  the  literature,  there  are  many  studies  of  battery  model¬ 
ing  [13,14,16,20-28].  Modeling  of  reaction  kinetics  and  transport 
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Nomenclature 

dm 

specific  surface  area  of  Fe  (cm-1 ) 

as 

specific  surface  area  of  FeCl2  (cm-1) 

Ca 

concentration  of  NaAlCl4  (mol cm-3) 

Cb 

concentration  of  NaCl  (mol  cm-3) 

^r,b 

bulk  concentration  of  ferrous  complex  (mol  cm-3) 

Cr,e 

equilibrium  concentration  of  ferrous  complex 
(mol  cm-3) 

*T,bsat 

saturation  bulk  concentration  of  ferrous  complex 
(mol  cm-3) 

Cr.esat 

saturation  equilibrium  concentration  of  ferrous 
complex  (mol  cm-3 ) 

Cr.s 

surface  concentration  of  ferrous  complex 
(mol  cm-3) 

CT 

total  concentration  (mol cm-3) 

D 

diffusion  coefficient  of  electrolyte  (cm2  s-1 ) 

De 

effective  diffusion  coefficient  of  electrolyte 
(cm2  s-1) 

F 

Faraday’s  constant  (C  mol-1 ) 

H 

height  of  the  cell  (cm) 

I 

apparent  current  density  at  separator  (A cm-2) 

*0 

exchange  current  density  (A cm-2) 

u 

superficial  current  density  in  matrix  phase  (A cm-2) 

h 

superficial  current  density  in  electrolyte  phase 

(A  cm-2) 

j 

local  transfer  current  (A cm-3) 

/Qp.FeCl 

mole  fraction  equilibrium  constant  for  the  solubility 
ofFeCh 

^sp.NaCl 

solubility  product  of  NaCl  (mol2  cm-6) 

/Co,  K2, 

K M  mole  fraction  equilibrium  constants  for 
AlCl3-NaCl  solvent  equilibrium 

km 

mass  transfer  coefficient  of  ferrous  complex 
between  Fe  and  bulk  (cm  s-1 ) 

ks 

mass  transfer  coefficient  of  ferrous  complex 
between  FeCl2  and  bulk  (cm  s-1 ) 

kp 

rate  constant  for  NaCl  precipitation/dissolution 
reaction  (cm3  mol-1  s-1 ) 

Nr 

flux  of  ferrous  complex  (mol  cm-2  s-1 ) 

R 

gas  constant  (J  mol-1  K-1 ) 

R¥eC\2  p 

precipitation/dissolution  rate  of  FeCl2 

(mol  cm-3  s-1) 

/^NaCl  p 

precipitation/dissolution  rate  of  NaCl 

(mol  cm-3  s-1) 

r 

radial  distance  from  the  center  of  current  collector 
(cm) 

to 

outer  radius  of  the  current  collector  (cm) 

Fa 

outer  radius  of  the  separator  (cm) 

tc 

outer  radius  of  the  negative  electrode  (cm) 

*L 

outer  radius  of  the  positive  electrode  (cm) 

ts 

outer  radius  of  the  electrolyte  reservoir  (cm) 

T 

temperature  (K) 

t 

time  (s) 

transference  number  of  A1C14-  relative  to  the  com¬ 
mon  ion  velocity 

*2 

transference  number  of  Cl-  relative  to  the  common 
ion  velocity 

t* 

transference  number  of  Na+  relative  to  the  molar- 
average  velocity 

VA 

molar  volume  of  molten  NaAlCl4  salt  (cm3  mol-1 ) 

VB 

molar  volume  of  molten  NaCl  salt  (cm3  mol-1 ) 

Ve 

molar  volume  of  electrolyte  (cm3  mol-1 ) 

I/Fe 

molar  volume  of  Fe  (cm3  mol-1 ) 

VFeci2  molar  volume  of  FeCl2  (cm3  mol-1 ) 

Wiaci  molar  volume  of  NaCl  precipitate  (cm3  mol-1 ) 

V  cell  potential  (V) 

V0c  open-circuit  cell  potential  (V) 

v  molar-average  electrolyte  velocity  (cm  s-1 ) 

xA  mole  fraction  of  NaAlCl4 

xAsat  saturation  mole  fraction  of  NaAlCl4 

Xb  mole  fraction  of  NaCl 

*aici3,  ^ai2ci6 *  *ai2ci-’  xFe(Aici4)2-  mole  fractions  of  A1C13, 
A12C16,  A12C17-  and  Fe(AlCl4)42- 

Greek  letters 

aa,  oi c  anodic  and  cathodic  transfer  coefficients 

yA  NaAlCl4  activity  coefficient 

£  porosity 

£pe  volume  fraction  of  Fe 

£FeC1  volume  fraction  of  FeCl2 

£NaCi  volume  fraction  of  NaCl  precipitate 

r]  total  overpotential  (V) 

k  Electrolyte  conductivity  (S  cm-1 ) 

Ke  Electrolyte  effective  conductivity  (S  cm-1 ) 

a  Iron  conductivity  (S  cm-1 ) 

ex e  Iron  effective  conductivity  (S  cm-1 ) 

0!  potential  in  matrix  phase  (V) 

(f>2  potential  in  electrolyte  phase  (V) 


processes  in  the  cell  is  critical  to  predict  the  change  of  cell  poten¬ 
tial  with  the  depth  of  discharge  and  to  evaluate  the  outcomes  of 
changes  in  design  parameters  [14,28].  Although  there  are  many 
studies  in  the  modeling  of  secondary  lithium  batteries  [20-24], 
the  modeling  of  sodium-metal  chloride  batteries  is  more  limited 
[13,14,16,28].  Sudoh  and  Newman  [14]  discuss  a  very  detailed 
model  of  a  discharge-charge  cycle  of  Na/p"-Al203/NaAlCl4/FeCl2 
battery  based  on  the  macroscopic  theory  of  porous  electrodes  [29] 
and  concentrated  solution  theory  [30].  In  their  model,  the  precip¬ 
itation/dissolution  rate  of  NaCl  is  taken  into  account.  In  addition, 
the  mass  transfer  of  the  soluble  ferrous  complex  is  included  in  the 
electrode  reaction  rate  [14],  however  they  do  not  allow  for  redis¬ 
tribution  within  the  cell  via  transport  of  the  iron  species  through 
the  electrolyte. 

Bloom  et  al.  [  1 6]  simulated  the  discharge  of  sodium-nickel  chlo¬ 
ride  cells  with  a  model  that  does  not  include  the  solubility  of  NiCl2 
and  NaCl.  Orchard  and  Weaving  [28]  also  published  a  model  on  the 
discharge  of  sodium-iron  chloride  cells.  In  this  study,  the  solubil¬ 
ity  of  FeCl2  and  NaCl  are  not  considered  [28].  Vallance  and  White 
[13]  modified  Sudoh  and  Newman’s  model  [14]  and  created  a  two- 
dimensional  model  for  a  fluted  (3"-alumina  tube. 

In  this  paper,  a  mathematical  model  of  the  porous  cathode  of 
a  Na/p"-Al203/NaAlCl4/FeCl2  battery  during  a  discharge-charge 
cycle  is  presented.  The  cathode  is  modeled  using  the  macroscopic 
theory  of  porous  electrodes  [29].  Transport  processes  are  modeled 
using  Pollard  and  Newman’s  [30]  concentrated-solution  theory  for 
a  mixture  of  two  binary  molten  salts  in  a  porous  electrode.  Although 
the  previous  models  are  successful  in  defining  the  kinetics  and 
transport  in  the  cathode,  none  of  them  can  predict  the  movement 
of  the  metal  that  takes  place  in  the  cell  with  increased  cycling 
[9,1 2,1 5,18,19].  An  advance  offered  by  this  model  is  that  it  accounts 
for  the  change  in  the  solubility  of  FeCl2  within  the  cell  and  predicts 
the  relocation  of  the  iron  by  modeling  the  transport  assuming  iron 
is  dilute. 
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Molten  Electrolyte  (Aid. ,  Cl  .  Na',  AI.CI,,  Aid.,  AI.CU .  Fe(AICI4|,z  ) 

©  /  o 


ro  Fe  FeCI2  NaCI  rL  rs  rA  rc 


Fig.  1.  Schematic  diagram  of  the  sodium-iron  chloride  battery  in  the  model. 


2.  Model  development 

The  isothermal,  constant  current  discharge-charge  of  a  Na-FeCl2 
battery  is  represented  using  a  one-dimensional,  cylindrical  model. 
As  seen  in  Fig.  1,  the  cell  is  composed  of  six  parts:  the  cathode 
current  collector  (r<r0),  the  positive  porous  Fe/FeCl2  electrode 
(between  r0  and  rL),  the  sodium  tetrachloroaluminate  reservoir 
(between  rL  and  rs),  the  p"-alumina  solid  electrolyte  (between  rs 
and  rA),  the  negative  liquid  sodium  electrode  (between  rA  and  rc) 
and  the  anode  current  collector  [14].  In  this  study,  only  the  positive 
porous  electrode  is  modeled. 

The  positive  porous  electrode  (between  r0  and  rL)  is  composed 
of  a  matrix,  which  consists  of  NaCI  crystals  and  porous  iron  particles 
partially  coated  with  FeCl2,  and  a  molten  electrolyte,  a  mixture  of 
AICI3  and  NaCI.  Transport  equations  are  derived  from  Pollard  and 
Newman’s  [30]  study  for  a  mixture  of  two  binary  molten  salts  with 
a  common  ion  in  a  porous  electrode.  As  in  their  study,  A  and  B  are 
used  for  NaAlCl4  and  NaCI  salts  and  1,  2  and  3  are  used  for  A1C14-, 
Cl-  and  Na+  ions,  respectively. 

The  electrolyte  is  a  concentrated  solution  of  A1C14-,  Cl-  and 
Na+  [30].  From  solvent-equilibria  studies  of  AlCl3-NaCl  melts,  it  is 
known  that  the  electrolyte  also  contains  Al2Cl6,  AICI3  and  Al2Cl7- 
ions  in  low  concentrations  [3 1  ].  In  addition  to  these  ions,  there  is  the 
soluble  ferrous  complex,  most  probably  in  the  form  of  Fe(AlCl4)42_ 
[14,17].  The  concentrations  of  this  ferrous  complex  on  FeCl2  and 
Fe  surfaces  and  in  the  bulk  may  not  be  the  same  depending  on  the 
interfacial  reaction  rates.  The  schematic  diagram  of  the  Fe/FeCl2 
electrode  defining  the  equilibrium,  bulk  and  surface  concentrations 
of  the  soluble  ferrous  complex  can  be  seen  in  Fig.  2. 


Fig.  3.  Solubility  of FeCh,  cr,e,  as  a  function  of  NaAlCl4  mole  fraction,  xA,  for  different 
values  of  I<s PjFeci  in  a  semi-log  plot. 


2.1.  Solubility  of  FeCl2 

The  equilibrium  concentration  of  the  complex  on  the  FeCl2  sur¬ 
face  is  dictated  by  the  solubility  of  FeCl2.  In  the  previous  models 
[13,14],  it  was  assumed  that  the  solubility  concentration  of  iron 
chloride,  cr,e,  is  constant  within  the  cell.  In  this  model,  we  allow  it 
to  change  within  the  cell  as  a  function  of  radial  position  and  time. 

In  a  previous  study  [28],  the  solvent  equilibrium  of  AlCl3-NaCl 
melts  is  described,  and  three  mole  fraction  equilibrium  constants, 
K 0,  I<2,  and  I<M  were  found.  In  this  study,  another  equilibrium  con¬ 
stant,  Ksp.FeCb  is  defined  for  the  solubility  of  FeCl2.  Therefore,  the 
equilibrium  is  now  defined  with  the  reactions  (2)— (5): 


2A1C13(1)^  A12C16(1)  K 0 

(2) 

AICI4  -  +  AICI3  ^  Al2  Cl7  -  k2 

(3) 

2aici4-+*  Ai2ci7-+cr  km 

(4) 

FeCl2(s)-+2Al2Cl7-4>  Fe(AlCl4)42-  Ksp.FeCl 

(5) 

The  solubility  of  FeCl2,  cr>e,  which  is  determined  by  these  four 
equilibrium  reactions,  is  only  a  function  of  the  NaAlCl4  mole  frac¬ 
tion,  xA,  and  the  assumed  solubility  constant,  /C^Feci- In  Fig-  3,  the 

solubility  as  a  function  of  xA  is  shown  for  Ksp.Feci  values  between 
104  and  108.  It  can  be  seen  that,  as  /Csp.Feci  increases,  the  solubil¬ 
ity  of  FeCl2  also  increases.  The  model  predicts  a  significant  change 
in  the  solubility,  especially  for  xA  values  between  0.6  and  0.9.  The 
assumed  relationships  for  cr,e  are  given  in  Appendix  A. 


Molten  Electrolyte 
(AICI4-,  Cl',  Na+,  AI2CI6, 
AICI3,  AI2CI/,  Fe(AICI4)42') 


2AICI3(I)  =  AI2CI6(I) 

AICI4-  +  AICI3  =  AI2CI7- 
2AICI4'  =  AI2CI7-  +  Cl 

2(s) +  2AI2CI7'  =  Fe(AICI4)42 


FeCL  Surface 


K, 


sp,FeCI 


Fe  Surface 


Fig.  2.  Schematic  diagram  of  the  Fe/FeCh  electrode  showing  the  equilibrium  (cr>e),  bulk  (cr>b)  and  surface  (cr>s)  concentrations  of  the  soluble  ferrous  complex,  Fe(AlCl4)42~. 
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2.2.  Concentration  of  the  soluble  ferrous  complex  in  the 
electrolyte 

The  redistribution  of  iron  in  the  cell  as  a  result  of  cycling  indi¬ 
cates  there  is  a  movement  of  the  soluble  ferrous  complex  within 
the  cell.  For  this  reason,  the  flux  of  Fe(AlCl4)42_  was  included  in 
the  model.  The  electrolyte  is  a  concentrated  solution  composed  of 
AlCLr,  Cl",  Na+,  A12C16,  A1C13,  A12C17-  and  Fe(AlCl4)42-  ions.  Since 
the  concentration  of  the  species  1,  2  and  3  are  very  high  compared 
to  the  other  ions,  in  the  derivation  of  the  flux  equation  for  the  fer¬ 
rous  complex,  the  presence  of  A12C16,  A1C13  and  AI2CI7-  ions  were 
neglected.  Thus  the  multicomponent  diffusion  equation  [32  ]  for  the 
ferrous  complex,  Eq.  (6),  only  contains  species  1,2,3  and  the  ferrous 
complex,  denoted  with  the  subscript  r: 


2.4.  Reduction! oxidation  reaction  rate 

The  reaction  that  takes  place  at  the  cathode  is  shown  in  Eq.  ( 1 3 ): 

FeCl2(S)  +  2e  -<-»■  Fe(S)  +  2Cl  (13) 

The  reaction  rate  accounts  for  the  mass  transfer  of  ferrous  com¬ 
plex  from  the  bulk  to  the  Fe  surface  as  described  in  Eq.  (14).  This 
rate  expression  is  a  modification  of  Eq.  (12)  in  Ref.  [14],  with  the 
appropriate  changes  for  estimation  of  cr,e  and  cr  b.  All  the  other  mass 
transfer  rates  and  areas  are  the  same  with  the  previous  studies 
[13,14], 

.  =  exp((o!gF/Rr)?7)  -  ((cr>b/cr,e)exp(-(qcF/RT)?7)) 

(1/io^m)  +  (l/2Fcrie)  (l/fcm^m)  QXp(—(otcF/RT)f]) 


Cr,bVM r=^  -  Vr)+C-^(V2-Vr)+C-^(V3-V^  (6) 

where  / xr  is  the  electrochemical  potential  of  the  ferrous  complex, 
cT  is  the  total  concentration,  crb,  c\ ,  c2  and  c3  are  the  concentrations 
of  the  species  and  vr,  iq,  v2  and  i/3  are  the  velocities  of  the  species. 

Since  we  do  not  have  independent  measurement  of  the  diffusion 
coefficients  of  the  species,  Dn,  we  assume  Dr  1  =Dr2  =Dr3  =De.  Also 
since  cr b  «  C\ ,  c2  and  c3,  we  assume: 

cT  ^  Cl  +  c2  +  c3  (7) 

and 


CjV*  ~  N-[  +  N2  +  N3  (8) 

where  v  is  the  molar-average  velocity  and  N\ ,  N2  and  N3  are  the 
fluxes  of  the  species. 

With  these  assumptions,  Eq.  (6)  is  rewritten  to  give  the  flux  of 
ferrous  complex,  Nr: 


Nr  =  Cr>bl/r  =  - 


De 

RT 


Cr,bV/Xr  +  Cr,bV* 


(9) 


Inserting  the  Gibbs-Duhem  equation,  using  the  definition  of 
the  effective  diffusion  coefficient  of  the  electrolyte,  De  =  Dc15,  and 
neglecting  electrical  migration  since  the  transference  number  of 
the  iron  species  is  essentially  zero,  Eq.  (10)  is  obtained: 


where  the  total  overpotential,  p,  is  given  by: 

4  =  01-02  (15) 

At  any  point  in  the  cell  where  the  volume  fraction  of  iron  chlo¬ 
ride  is  zero  during  the  discharge,  the  transfer  current  is  set  to  zero. 


2.5.  Precipitation  rate  ofNaCl 


The  second  reaction  taking  place  inside  the  porous  cathode  is 
the  precipitation/dissolution  reaction  of  NaCl,  which  is  given  by 
Eq.  (16).  The  rate  of  this  reaction  is  shown  in  Eq.  ( 1 7)  [  1 4].  When  Xa 
is  equal  to  its  saturation  value,  the  rate  is  equal  to  zero.  When  it  is 
lower  than  the  saturation  value,  RNaC1  p  is  positive  indicating  there 
is  precipitation  of  NaCl  in  the  cell. 


Na+  +  Cl-  **  NaCl(s) 

^NaCl  p  =  kp  ^  ^ 2  l^sp.NaCl^ 

where  the  average  molar  volume  of  the  electrolyte,  Ve: 
Ve=(VA-VB)xA  +  VB 


(16) 

(17) 


(18) 


2.6.  Precipitation  rate  ofFeCh 


Nr  =  — D£15Vcr  b  +  crbv*  (10) 

A  material  balance  of  the  soluble  ferrous  complex  in  the 
electrolyte  enables  calculation  of  the  bulk  concentration  of 
Fe(AlCl4)42-,  crb.  When  the  quasi  steady-state  assumption  is 
applied,  the  mass  transfer  rate  of  ferrous  complex  from  the  FeCl2 
surface  to  the  bulk  is  equal  to  the  flux  of  ferrous  complex  in  the  elec¬ 
trolyte  and  the  electrochemical  reaction  rate.  The  material  balance 
for  crb  is  shown  in  Eq.  (1 1 ): 


The  precipitation  rate  of  iron  chloride  is  given  by: 

^FeCl2p  =  -Ms^e  -  Cr?b)  (19) 

The  precipitation/dissolution  rate  is  zero  when  the  equilibrium 
and  bulk  concentrations  of  the  ferrous  complex  are  equal  to  each 
other.  The  rate  is  positive  showing  there  is  precipitation  of  FeCl2 
when  the  bulk  concentration  of  the  ferrous  complex  is  higher  than 
its  equilibrium  concentration. 


0  —  —  V  •  Nr  +  —  +  /<s<3s(Cr,e  _  *T,b)  (1 1 ) 

where  j  is  the  local  transfer  current,  ks  is  the  mass  transfer  coef¬ 
ficient  of  ferrous  complex  between  FeCl2  and  bulk  and  as  is  the 
specific  surface  area  of  FeCl2. 

2.3.  Concentration  of  the  soluble  ferrous  complex  on  the  Fe 
surface 

The  surface  concentration,  cr,s,  is  calculated  by  equating  the 
mass  transfer  rate  of  ferrous  complex  from  the  bulk  to  the  Fe  surface 
to  the  electrochemical  reaction  rate  as  given  in  Eq.  (12): 

=  —  kmQm(Cr,b  —  cr,s)  (12) 

where  km  is  the  mass  transfer  coefficient  of  ferrous  complex 
between  Fe  and  bulk  and  am  is  the  specific  surface  area  of  Fe. 


2.7.  Material  balances  on  iron ,  iron  chloride  and  sodium  chloride 


The  equations  used  for  the  calculation  of  the  volume  fractions  of 
iron,  iron  chloride  and  sodium  chloride  are  shown  in  Eqs.  (20),  (21 ) 
and  (22),  respectively.  The  total  solids  porosity,  s,  given  in  Eq.  (23) 
is  calculated  based  on  the  fact  that  the  summation  of  the  porosity 
and  the  volume  fraction  of  the  matrix  is  equal  to  1. 


d£pe 

dt 

d£FeCi2 

~~dF~ 

d£NaC\ 

dt 


=  WeCl2^FeCl2  p 
=  ^NaCl^NaCl  p 


(20) 

(21) 

(22) 


9t=2Fj~W<P 


^,2  A  -  ^sp.NaCl  1  +WeCI2^sas(cr,e  —  cr,b)  (23) 
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2.8.  Material  balance  on  electrolyte 


2.13.  Cell  potential 


The  mole  fraction  of  NaAlCU,  xA,  is  the  variable  in  the  model 
used  to  set  the  composition  of  the  electrolyte.  A  material  balance 
is  the  same  as  given  in  Ref.  [13]: 

=  ^eXA,?NaC|p  +  ^eXAp  ~ v*VXa  +  v ' (°£l'5 Vxa) 

-  D£a-5  Va  T  Vb  (Vxa )2  +  ^  VxA  (24) 

ve  2t 


2.9.  Current  densities  in  electrolyte  and  matrix  phases 


One  of  the  most  important  features  of  battery  modeling  is  the 
prediction  of  the  external  cell  potential.  Ref.  [14]  discusses  the  cal¬ 
culation  of  the  terminal  voltage  in  detail.  In  this  paper,  a  simplified 
version  of  their  equation  is  used  under  the  assumptions  of  constant 
reservoir  and  ceramic  electrolyte  resistances  and  constant  nega¬ 
tive  electrode  overpotential.  With  these  assumptions,  the  change 
in  potential  with  time  can  simply  be  calculated  using: 

V  =  V0c  +  (</-i)r=ro-(</>2)r=rL+(^lnT^^)r  ^  (31) 

where  V0c  is  the  open-circuit  cell  potential  and  0i  and  02  are  the 
potentials  in  the  matrix  and  electrolyte  phases,  respectively. 


Ohm’s  law  is  used  to  define  the  current  densities  in  the  matrix 
and  electrolyte  phases,  i\  and  i2 ,  that  are  shown  in  Eqs.  (25)  and 
(26),  respectively  [  1 3,1 4,30].  The  summation  of  these  two  variables 
is  equal  to  the  apparent  current  density,  I  (measured  at  rs)  at  any 
point  in  the  cathode  [14].  Effective  conductivities  of  the  metal  and 
electrolyte  are  defined  as  ae  =  as 15  and  /ce  =  [13,14].  In  the 

derivations,  the  transference  number  of  sodium  ion,  1 1  is  assumed 
to  be  0.5  since  the  sodium-ion  concentration  is  half  of  the  total 
concentration,  and  t\  and  tS,  are  taken  as  xA  and  Xb,  respectively 
[14]. 


2.14.  Initial  conditions 

Since  the  cell  is  fully  charged  before  the  first  discharge,  the  ini¬ 
tial  conditions  for  £,  £Fe,  £Feci  and  %aci*  are  0-546,  0.184,  0.26  and 
0.01,  respectively  [13,14].  For  the  mole  fraction  of  A,  xA,  the  satu¬ 
ration  value  of  0.8972  [13]  is  used.  It  is  assumed  that  j  is  constant 
throughout  the  cell  initially,  and  the  initial  conditions  for  the  other 
5  variables,  <P\,  @2,  h,  h  and  p,  are  calculated  using  this  assump¬ 
tion.  Finally,  Eq.  (32),  which  is  derived  using  the  fluxes  of  species  1, 
2  and  3,  is  used  for  the  initial  condition  of  v  [14]. 


i\  =  — cre  V0t 


(25) 


h 


fVtl+^+  RTt1 

ae+Ke{  <7er  F(  l-xA)xA 


rflnyA 

dlnxA 


(26) 


,  Vre-2VrB-VFe  +  VrFea 

V  =  - -In 

2  F  2 

2.15.  Boundary  conditions 


(32) 


h+h  =  Y  (27) 

where  0i  is  the  potential  in  the  matrix  phase  and  yA  is  the  activity 
coefficient. 

2. 1 0.  Local  transfer  current 

The  definition  of  the  local  transfer  current  is  [13,14]: 
j  =  V  •  i2  (28) 

2.11.  Molar  average  velocity 

The  change  in  the  porosity  with  time  and  position  creates  a 
velocity  field  inside  the  cell,  which  is  given  in  Eq.  (29)  [13,14].  This 
velocity  within  the  cell  results  in  the  convective  mass  transfer  of 
the  species,  therefore  it  must  be  considered  in  the  model. 

v  ■  v*  =  -VFe^F2VBj  +  (V^NaCl  -  Vb)  %>Clp 
+  (Va  -  VB)  V  ■  [Ds15(ca  +  cB)VxA] 

+  ■  (t;;2)-  ^ V  ■  (tp2)  -  VFea2ksas(Cr,e  -  cr, b)  (29) 

2.12.  Total  iron  amount 

In  order  to  describe  the  movement  of  the  iron  in  the  cell,  a 
parameter  for  the  total  iron  amount  is  used  to  show  iron  redis¬ 
tribution  within  the  battery: 

Total|ron  =  p-  +  +  Crb  (30) 

^Fe  ^FeC^ 


The  boundary  conditions  at  r0  are  given  in  Eqs.  (33)-(39): 


i2  =  0 

(33) 

-r 

ro 

(34) 

V &2  =  0 

(35) 

_  rsI 

V  p  = - — 

aer0 

(36) 

< 

> 

II 

0 

(37) 

v*  =  0 

(38) 

Vcrb  =  0 

(39) 

The  boundary  conditions  at  rL  are: 

i  /rs 

l2  =  TL 

(40) 

h  =  0 

(41) 

01  =0 

(42) 

0 

II 

< 

> 

(43) 

,  Ve-2VB-VFe  +  VFeci2. 

V  =  2  F  '2 

(44) 

O 

II 

XJ 

u 

> 

(45) 

The  14  variables  in  the  model,  <2>i,  @2 ,  i'i,  i2,j ,  p,  xA,  £,  £Fe,  £Feci* 
£Naci»  v*y  cr,e  and  cr  b,  are  calculated  solving  the  Eqs.  (11),  (14),  (15), 
(20)-(29)  and  (A6)  using  a  block  tri-diagonal  matrix  algorithm  in 
FORTRAN  [32].  The  simulations  were  carried  out  with  501  node 
points  and  time-step  sizes  of  20.61  s  during  discharge  and  9  s  dur¬ 
ing  charge.  For  continuous  cycling,  time-step  sizes  of  6  s  and  9  s  are 
used  for  discharge  and  charge,  respectively.  In  order  to  test  the  con¬ 
vergence,  time-step  size  was  halved,  resulting  in  negligible  changes 
in  the  results.  In  a  similar  manner,  doubling  of  node  points  did  not 
change  the  results  to  any  appreciable  extent. 
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Table  1 

Parameters  in  the  model. 


r0 

0.25  cm 

rL 

2.5  cm 

^Middle 

1.3705  cm 

rs 

2.8  cm 

H 

30  cm 

T 

573  K 

I 

-30  mA  cm-2  discharge/10  mA  cm-2  charge 

VFe 

7.1  cm3  mol-1 

YFeC12 

40.1  cm3  mol-1 

^NaCl 

27.0  cm3  mol-1 

VA 

121.6  cm3  mol-1 

Vb 

37.06  cm3  mol-1 

kp 

0.1  cm3  mol-1  s_1 

/Csp.NaCl 

8.06  x  10-6  mol2  cm-6 

aa=ac 

1 

D 

5.135  x  10“6  cm2  s-1 

cr 

3.5  x  104  Son-1 

/Co 

3.85  xlO5 

k2 

103 

I<m 

2.37  x  10“6 

Vqc 

2.32  V 
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3.  Results  and  discussion 


Fig.  5.  The  change  of  NaAlCU  mole  fraction,  xA,  with  radial  position  at  different 
DODs  during  discharge  (/CsPiFe a  =  106). 


Results  are  discussed  assuming  an  operating  temperature  of 
300  °C.  The  values  of  parameters  used  in  the  simulations  are  shown 
in  Table  1.  Most  of  these  values  are  taken  from  previous  work  for 
sodium-iron  chloride  battery  modeling  [13,14,31  ].  For  the  electri¬ 
cal  conductivity  and  activity  coefficient  equations,  the  equations 
in  the  Appendix  of  Ref.  [14]  were  used.  The  full  discharge  time  is 
calculated  as  46,040  s  (12.8  h)  with  -30  mAcnrr2  constant  current 
discharging  and  the  depth  of  discharge  (DOD)  is  computed  as  the 
ratio  of  the  actual  discharge  time  to  the  full  discharge  time.  The 
charge  time  is  taken  as  1 8,000  s  (5  h)  in  the  simulations.  In  all  of  the 
results  given  for  the  model,  I<sp  Feci  is  taken  as  1 06,  unless  otherwise 
stated. 

3.1.  Solubility  ofFeCl2  and  bulk  concentration  of  the  ferrous 
complex 

The  change  in  iron  chloride  volume  fraction,  NaAlCl4  mole  frac¬ 
tion  and  solubility  of  FeCl2  within  the  cell  during  discharge  can  be 
seen  in  Figs.  4-6,  respectively.  Initially  (DOD  =  0),  the  solubility  is 


constant  along  the  cell  since  xA  is  equal  to  its  saturation  value  at 
every  point  in  the  cell.  With  increasing  time  (DOD  =  0.2),  the  sol¬ 
ubility  begins  to  decrease  at  the  electrode-reservoir  interface,  rL. 
This  is  expected  since  xA  has  the  same  trend  at  DOD  =  0.2  (Fig.  5) 
with  the  minimum  value  at  rL.  For  larger  discharge  times  (DOD  =  0.7 
and  0.9),  the  minimum  value  of  xA  shifts  inward  to  r0  (Fig.  5)  and 
the  solubility  decreases  throughout  the  cell  with  the  same  trend 
(Fig.  6). 

The  bulk  concentration  of  the  ferrous  complex  with  respect 
to  position  and  time  is  shown  in  Fig.  7.  As  it  can  be  seen  in  the 
figure,  at  low  discharge  times  (DOD  =  0.2)  the  bulk  concentration 
follows  the  equilibrium  concentration;  it  decreases  with  increasing 
r.  As  discharge  time  increases  (DOD  =  0.7),  although  iron  chloride  is 
depleted  near  the  electrode-reservoir  interface  (Fig.  4),  the  con¬ 
centration  of  the  soluble  ferrous  complex  in  the  electrolyte  is  not 
zero.  This  result  shows  that  the  flux  of  the  ferrous  complex  becomes 
more  significant  with  increasing  depth  of  discharge  in  part  because 
FeCl2(S)  is  no  longer  present  to  buffer  variations  in  the  bulk  iron 
concentration.  As  the  discharge  time  increases  further  (DOD  =  0.9), 


Fig.  4.  The  change  of  iron  chloride  volume  fraction,  eFea>  with  radial  position  at 
different  DODs  during  discharge  (/CspFeci  =  106). 


Fig.  6.  The  change  of  solubility  of  FeCl2,  cr,e,  with  radial  position  at  different  DODs 
during  discharge  (/CSp,Feci  =  106). 


D.  Eroglu,  AC.  West  /  Journal  of  Power  Sources  203  (2012)211  -221 


217 


Fig.  7.  The  change  of  bulk  concentration  of  ferrous  complex,  cr  b,  with  radial  position 
at  different  DODs  during  discharge  (I<s PiFeci  =  106). 


the  electrolyte  concentration  of  iron  approaches  zero  near  the  rL 
boundary.  At  this  high  depth  of  discharge,  iron  chloride  is  depleted 
in  the  majority  of  the  cell  (Fig.  4),  and  the  electrolyte  concentration 
of  the  ferrous  complex  is  very  low,  even  near  r0.  Therefore,  the  flux 
of  the  ferrous  complex  is  negligibly  small;  it  is  not  enough  to  create 
non-zero  bulk  concentrations  near  rL.  The  discontinuities  seen  in 
Fig.  7  (DOD  =  0.7  and  0.9)  occur  at  the  points  where  the  iron  chlo¬ 
ride  phase  is  calculated  to  disappear.  We  have  confirmed  that  the 
discontinuities  do  not  propagate  numerical  errors. 

3.2.  Relocation  of  iron  within  the  cell 

The  change  in  the  total  iron  amount  at  r0,  rMiddle  and  rL  during  a 
discharge-charge  cycle  is  shown  in  Fig.  8.  In  the  figure,  it  can  be  seen 
that  there  is  a  sudden  increase  in  the  total  iron  amount  at  rL  dur¬ 
ing  the  early  stages  of  discharge.  This  increase  in  total  iron  amount 
can  be  explained  by  a  sharp  increase  in  Fe  concentration.  This  sharp 
increase  shows  that  there  is  a  significant  flux  of  the  ferrous  complex 
present  towards  rL  during  the  early  stages  of  the  discharge.  After 
this  initial  increase,  the  total  iron  content  does  not  change  at  rL  until 
the  end  of  discharge.  As  the  iron  chloride  is  depleted  near  rL,  the 
reaction  front  migrates  inward  towards  r0,  and  the  increase  in  Fe 
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Fig.  8.  The  change  of  total  iron  amount  with  time  at  different  radial  positions 

(Ksp, Fed  =  106  )• 
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Fig.  9.  Comparison  of  bulk  concentration  results  for  the  complete  model  (model 
4  -  flux,  variable  solubility)  and  different  simplifications  (model  1  -  no  flux,  con¬ 
stant  solubility),  (model  2  -  flux,  constant  solubility),  (model  3  -  no  flux,  variable 
solubility).  For  all  cases  DOD  =  0.7  and  I<sp, Feci  =  iO6. 


concentration  also  ends.  This  suggests  that  unlike  the  initial  stages, 
the  flux  of  the  ferrous  complex  is  not  enough  to  create  a  change 
in  Fe  concentration  or  total  iron  amount.  As  discharge  proceeds, 
cr  b  also  goes  to  zero  and  the  total  iron  amount  reflects  only  the 
Fe  concentration  at  rL.  At  rMiddle  and  r0,  the  iron  amount  decreases 
slightly  during  discharge.  These  results  show  that  there  is  a  net 
flux  of  ferrous  complex  from  r0  to  rL.  The  flux  of  the  ferrous  com¬ 
plex  is  reversed  during  charging  of  the  cell,  causing  the  total  iron 
amount  to  decrease  at  rL  and  increase  at  rMiddle  and  r0.  As  a  result 
of  a  discharge-charge  cycle,  the  total  amount  of  iron  is  reduced 
slightly  at  rL.  This  small  change  in  the  total  iron  amount  as  a  result 
of  one  cycle  may  become  significant  after  several  dozen  cycles. 


3.3.  Significance  of  solubility  ofFeCl2  variation  and  flux  of  the 
ferrous  complex  in  the  model 

In  order  to  examine  the  role  of  the  iron  chloride  solubility  change 
and  ferrous  complex  flux  in  the  prediction  of  iron  redistribution 
inside  the  cell,  four  model  predictions  are  compared.  In  the  first 
model,  the  solubility  of  iron  chloride  is  constant  throughout  the 
cell  at  its  saturation  value.  In  addition,  the  flux  of  the  ferrous  com¬ 
plex  is  not  included  in  the  material  balance  of  the  ferrous  complex 
in  the  bulk.  This  first  model  is  essentially  identical  to  Sudoh  and 
Newman’s  model  [13,14].  In  the  second  model,  the  solubility  of  iron 
chloride  is  constant  throughout  the  cell,  but  the  flux  of  the  ferrous 
complex  is  included.  In  the  third  model,  the  variation  in  FeCl2  solu¬ 
bility  is  taken  into  account,  but  the  flux  is  not  included.  Finally,  the 
fourth  model  is  the  one  described  in  this  paper;  counting  for  both 
the  solubility  change  and  ferrous  complex  flux  inside  the  cell.  In 
order  to  compare  these  four  models,  the  bulk  concentration  profile 
at  DOD  =  0.7  is  chosen  as  shown  in  Fig.  9. 

When  the  models  with  no  flux  (models  1  and  3)  are  compared 
with  the  models  accounting  for  the  flux  of  the  ferrous  complex 
(models  2  and  4)  in  Fig.  9,  it  can  be  seen  that  the  presence  of  the  flux 
in  the  model  resulted  in  prediction  of  higher  bulk  concentrations 
near  the  rL  boundary.  For  instance,  the  bulk  concentration  is  non¬ 
zero  for  model  4  near  the  ri  boundary  although  it  is  zero  for  model 
3.  From  these  results  it  can  be  concluded  that  the  effect  of  the  flux  of 
the  ferrous  complex  is  significant  at  moderate  DODs  near  rL,  where 
iron  chloride  is  depleted. 
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Fig.  11.  Comparison  of  iron  chloride  volume  fraction  profiles  for  different  /(sp,Feci 
values  during  discharge. 


Fig.  10.  Comparison  of  total  iron  amount  results  for  the  complete  model  (model 
4  -  flux,  variable  solubility)  and  different  simplifications  (model  1  -  no  flux,  con¬ 
stant  solubility),  (model  2  -  flux,  constant  solubility),  (model  3  -  no  flux,  variable 
solubility).  For  all  cases  DOD  =  0.7  and  I<sPt FeCi  =  106. 

The  effect  of  variation  of  FeCl2  solubility  on  the  electrolyte  con¬ 
centration  is  also  apparent  in  Fig.  9,  especially  near  r0,  where  solid 
iron  chloride  is  still  present.  Electrolyte  concentration  does  not 
change  near  r0  for  the  models  with  constant  solubility  (models  1 
and  2)  whereas  it  decreases  with  radial  distance  for  models  with 
variable  solubility  (models  3  and  4).  The  variation  in  cr,e  results  in 
lower  electrolyte  concentrations  near  this  boundary. 

In  Fig.  10,  the  total  iron  amount  at  rL  as  a  function  of  discharge 
time  for  these  four  models  are  compared.  This  figure  clearly  dis¬ 
plays  the  importance  of  the  flux  of  the  ferrous  complex  in  the 
redistribution  of  the  iron  inside  the  cell  since  models  2  and  4  exhibit 
a  significant  change  in  the  total  iron  amount.  The  increase  in  total 
iron  amount  predicted  by  model  2  is  higher  than  model  4  since  in 
model  2,  solubility  of  FeCl2  is  constant  at  its  maximum  value. 

It  can  be  concluded  that  although  the  variation  in  the  FeCl2  sol¬ 
ubility  is  important,  the  addition  of  the  flux  of  the  ferrous  complex 
into  the  model  has  a  more  important  effect  on  the  results. 

3.4.  Effect  ofKsp  FeCi  on  the  discharge-charge  cycle  of  the  cell 

The  solubility  product,  7Csp  Feci,  has  not  been  previously  reported 
in  the  literature.  Sudoh  and  Newman  [14]  assumed  a  constant  equi¬ 
librium  ferrous  complex  concentration  of  4.1  x  10-8molcm-3  in 
their  study.  This  value  of  cr>e  corresponds  to  Ksp.Feci  of  2.7  x  104 
calculated  at  the  saturation  value  of  xA  according  to  our  model. 
We  have  carried  out  simulations  for  assumed  values  of  104,  105, 
106, 107  and  108.  When  I<sp,Fec\  <  105,  model  predictions  are  not  in 
accord  with  the  previous  results  reported,  most  probably  because 
cFie  is  too  low  in  the  system,  seriously  affecting  the  kinetics.  The 
results  show  consistent  behavior  for  all  TCsp.Feci  values  only  when 
it  is  equal  or  higher  than  105.  I<sp >Feci  values  higher  than  108  are 
not  considered  since  the  main  assumption  in  our  model  is  that  iron 
chloride  is  sparingly  soluble  in  the  electrolyte  and  this  range  leads 
to  too  high  of  a  concentration. 

3.4A.  Effect  ofI<sp  FeCi  on  the  iron  chloride  volume  fraction 

Iron  chloride  volume  fractions  with  respect  to  position  and 
time  for  /(sp,Feci  values  of  105,  106,  107  and  108  are  shown  in 
Figs.  11  and  12,  for  discharge  and  charge,  respectively.  During  dis¬ 
charge,  for  all  TCsp.Feci  values,  the  same  trend  is  seen;  £Feci  decreases 
with  time  throughout  the  battery  and  it  becomes  depleted  for 


positions  approaching  rL  at  high  DODs  since  the  reaction  front 
migrates  inwards  towards  the  positive  current  collector  [12]. 
When  the  curves  for  different  solubility  products  are  compared, 
it  can  be  seen  that  there  is  a  significant  difference  in  the  results 
only  when  KspFeC\  is  equal  to  105  (Fig.  11).  The  results  for  the 
other  three  I<spFe a  values  are  very  similar.  The  difference  between 
the  105  curve  and  the  others  increases  with  increasing  DOD.  For 
/Csp.Feci  values  of  106  and  higher,  iron  chloride  is  depleted  faster 
than  KSp,Feci =  1  °5  near  the  rL  boundary  due  to  higher  reaction 
rates.  When  comparing  curves  with  the  same  DOD,  a  higher  local 
transfer  current  at  r¥  must  result  in  a  lower  local  transfer  current 
at  r0.  Therefore  the  iron  chloride  volume  fraction  is  lower  at  r0  for 
/Csp.Feci  =  105  relative  to  the  higher  TQp.Feci  values. 

In  contrast,  during  charge,  the  value  of  /CsPiFeci  appears  to  affect 
behavior  near  the  rL  boundary.  As  it  can  be  seen  in  Fig.  1 2,  there  is  a 
sudden  decrease  for  £Feci  at  rL  for  108.  This  result  suggests  that  the 
equilibrium  value  of  the  ferrous  complex  has  a  significant  effect  on 
£Feci  at  the  rL  boundary  when  I<sp  Feci  is  higher  than  107.  This  effect 
is  not  seen  during  discharge  at  the  rL  boundary  because  FeCl2(S)  is 
depleted  for  any  value  of  the  solubility  product. 

3.4.2.  Effect  ofI<spFea  on  the  cell  potential 

The  change  in  the  cell  potential  with  time  during  discharge  for 
different  values  of  KSpyec\  can  be  seen  in  Fig.  13.  The  trend  for  all 


Fig.  12.  Comparison  of  iron  chloride  volume  fraction  profiles  for  different  KSp,Feci 
values  during  charge. 
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Fig.  13.  Comparison  of  cell  potentials  for  different  I<sp,Fea  values  during  discharge. 

of  the  curves  is  the  same;  it  starts  with  a  slight  decrease  in  the 
cell  potential  at  small  DODs,  then  a  significant  continuous  decrease 
is  observed,  and  finally  at  very  high  DODs,  a  very  steep  decrease 
occurs.  All  four  curves  are  the  same  until  35,000  s.  After  this  point 
there  is  a  clear  difference  among  the  curves;  they  have  the  steep 
potential  decrease  at  different  times.  As  for  the  previous  results, 
^sp.Feci =  1°5  has  the  most  significant  distinction  from  the  others. 
The  reason  why  the  steep  potential  decrease  is  seen  earlier  for  lower 
Ksp  feci  values  can  be  explained  by  the  slower  kinetics  at  the  reac¬ 
tion  front  in  these  systems  due  to  the  low  electrolyte  concentration 
of  ferrous  complex.  The  reaction  rate  at  the  reaction  front  is  signif¬ 
icantly  lower  for  lower  I< sPiFeci  values,  especially  for  Ksp  FeCi  =  105. 
The  effect  of  I<sp  Feci  on  the  cell  potential  is  observed  only  at  very 
high  DODs,  suggesting  that  the  sharp  potential  decrease  occurs 
when  reaction  rates  at  the  reaction  front  decreases  to  very  low  val¬ 
ues.  It  can  be  seen  that,  in  terms  of  the  cell  potential,  increasing  the 
solubility  product  higher  than  107  does  not  have  a  great  impact  on 
the  simulated  cell  potential  except  near  the  end  of  discharge. 

3.5.  Effect  of  continuous  cycling 

Continuous  cycling  simulations  were  conducted  to  study 
whether  iron  redistribution  continued  beyond  the  first  cycle. 
Here,  the  cycling  conditions  are  taken  as  discharge  for  12,000  s  at 


Fig.  14.  The  change  of  discharge  cell  potential  with  time  for  each  cycle  during 
continuous  cycling  (KsPiFea  =  106). 


-30  mA  cm-2  and  charge  for  18,000  s  at  10  mA  cm-2.  The  effect  of 
continuous  cycling  was  investigated  for  five  cycles. 

In  Fig.  14,  the  change  of  cell  potential  with  time  during  dis¬ 
charge  is  given  for  each  cycle.  It  can  be  seen  that  under  these  cycling 
conditions,  the  cell  potential  is  predicted  to  decrease  even  after  5 
cycles.  In  addition,  with  cycling  the  trend  of  the  potential  curve  also 
changes;  there  is  a  steep  decrease  in  the  potential  at  first  followed 
by  a  slighter  continuous  decrease.  This  sudden  potential  decrease 
starts  to  occur  around  the  same  DOD  at  each  cycle.  As  a  result  of  the 
simulation,  it  has  been  calculated  that  after  the  5th  discharge,  the 
cell  potential  decreases  by  nearly  10%.  The  change  in  the  discharge 
potential  curve  and  decline  in  discharge  potential  with  continu¬ 
ous  cycling  has  been  reported  previously  [19].  The  redistribution 
of  iron  inside  the  cathode  with  continuous  cycling  plays  an  impor¬ 
tant  role  in  this  potential  loss.  However  it  should  be  kept  in  mind 
that  the  discharge/charge  current  densities  and  simulation  times 
also  have  a  great  impact  on  the  calculated  potential  loss.  Therefore, 
this  extreme  of  potential  decrease  may  not  be  seen  with  different 
cycling  conditions.  In  addition,  it  has  seen  that  cell  design  has  a  crit¬ 
ical  effect  on  the  calculated  potential  loss.  For  instance,  simulations 
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Fig.  15.  The  change  of  total  iron  amount  with  time  at  different  radial  positions  during  continuous  cycling  (Ks p  Feci  =  106). 
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suggest  that  the  initial  volume  fraction  ratio  of  iron  to  iron-chloride 
creates  a  significant  difference  in  the  calculated  potential  loss. 

The  change  in  total  iron  amount  at  r0,  rMiddle  and  rL  with  con¬ 
tinuous  cycling  can  be  seen  in  Fig.  15.  In  the  previous  section,  a 
single  discharge-charge  cycle  was  investigated  and,  as  a  result  of 
the  simulations,  it  was  seen  that  iron  amount  increases  at  rL  dur¬ 
ing  discharge  and  decreases  during  charge.  The  same  trends  can  be 
seen  for  each  individual  cycle  in  Fig.  1 5.  As  explained  in  the  previous 
section,  the  initial  increase  is  mainly  due  to  the  sharp  increase  in 
iron  volume  fraction  because  of  the  significant  flux  of  the  ferrous 
complex.  After  iron  chloride  is  depleted,  the  increase  in  the  total 
iron  amount  ends  and  the  curve  stabilizes.  During  charging,  the  fer¬ 
rous  complex  diffuses  away  from  rL,  causing  the  Fe  concentration  to 
decrease  more  than  the  FeCl2  amount  and  the  bulk  concentration 
to  increase.  The  figure  also  shows  that  with  increased  cycling,  the 
iron  depletion  at  rL  at  the  end  of  charging  becomes  more  and  more 
significant.  For  instance,  as  a  result  of  the  simulation,  the  total  iron 
amount  decreases  by  ~1%  at  rL  at  the  end  of  the  fifth  cycle.  This 
result  suggests  that  iron  may  be  depleted  at  rL  around  the  end  of 
500th  cycle  causing  perhaps  the  failure  of  the  cell.  The  change  in 
the  iron  amount  at  r0  and  rMiddle  with  continuous  cycling  is  not  as 
significant  as  the  change  in  rL. 

Since  the  molar  average  velocity  at  rL  is  not  zero,  there  is  a  net 
flux  of  the  soluble  ferrous  complex  into  the  sodium  tetrachloroalu- 
minate  reservoir  as  a  result  of  continuous  cycling.  Flowever  it  was 
calculated  that  the  amount  of  iron  that  is  lost  into  the  reservoir 
due  to  this  flux  is  negligible  compared  to  the  total  decrease  in  the 
iron  amount  at  this  boundary.  Therefore,  the  redistribution  of  iron 
within  the  electrode  must  be  the  main  reason  of  the  iron  loss  at  rL. 


depleted  near  r l.  When  the  relocation  of  iron  is  considered,  it  was 
concluded  that  the  net  movement  of  iron  is  from  r0  to  rL  during 
discharge,  whereas  the  reverse  happens  during  charge.  The  effect 
of  solubility  constant  /CsPiFeci  was  also  studied,  and  effects  increase 
with  increasing  solubility,  until  1 06.  Finally,  the  effect  of  continuous 
cycling  was  examined.  It  was  predicted  that  there  is  a  deficiency  of 
iron  at  rL.  As  a  result  of  the  5th  cycle,  there  is  nearly  a  1%  decrease 
in  the  total  iron  amount  at  the  electrode-reservoir  boundary. 

Appendix  A.  Calculation  of  the  equilibrium  concentration 
of  the  ferrous  complex 


The  equilibrium  concentration  of  the  ferrous  complex  is  calcu¬ 
lated  using  Eqs.  (A1)-(A6).  The  concentrations  of  A1C14_,  Cl-  and 
Na+  are  all  expressed  in  terms  of  xA. 


/CMX  a 

aic,3  -  K2(  1  -xA) 

(Al) 

xA12C16  =  ^0XAICI3 

(A2) 

Kmx2a 

Xa‘2C17  “  CT(1  -xA)Ve 

(A3) 

XFe(AlCl4)2-  =  ^sp,FeCl*Al2C1- 

(A4) 

2 

Cj  =  _ 

Ve(l  -*A12C16  -*A1C13  -xA12C1-  -XFe(AlCl4)2-^ 

(A5) 

Cr.e  =  cTXFe(AlCl4)2- 

(A6) 

3.6.  Summary  and  future  work 

Minimizing  the  transport  and  redistribution  of  iron  inside  the 
cathode  should  be  considered  in  the  cell  design  in  order  to  prevent 
the  power  loss  seen  as  a  result  of  continuous  cycling.  The  model 
proposed  in  this  study  is  able  to  predict  the  iron  redistribution 
inside  the  cell  as  a  function  of  design  parameters.  For  instance,  it 
has  been  seen  that  one  of  the  design  parameters  in  Zebra  cells,  the 
initial  ratio  of  iron  to  iron-chloride  volume  fraction,  has  a  significant 
effect  on  iron  redistribution  inside  the  cell.  Simulations  suggest  that 
increasing  the  initial  ratio  of  iron  to  iron-chloride  volume  fraction 
decreases  the  amount  of  iron  depleted  at  rL. 

The  results  reported  in  this  study  show  that  incorporation  of 
the  variation  of  iron  chloride  solubility  and  the  flux  of  the  fer¬ 
rous  complex  into  the  model  of  the  porous  cathode  is  important 
to  predict  the  redistribution  of  the  metal  inside  the  cell.  For  future 
work,  the  iron  chloride  solubility  as  a  function  of  xA  would  be  valu¬ 
able  to  measure  directly.  Furthermore  the  metal  distribution  inside 
the  cathode  at  different  positions  at  different  DODs  would  allow 
for  direct  testing  of  model  predictions.  In  addition,  extending  the 
model  proposed  in  this  paper  to  a  Na/NiCl2  cell  would  be  of  value 
since  it  is  known  that  Na/NiCl2  batteries  have  some  superior  prop¬ 
erties  over  Na/FeCl2  batteries,  such  as  higher  open-circuit  potential 
[4,12]. 

4.  Conclusions 

In  this  study,  a  mathematical  model  for  the  positive  cathode  in 
a  sodium-iron  chloride  cell  with  p"-alumina  ceramic  and  molten 
NaAlCU/NaCl  electrolyte  was  extended  by  accounting  for  variable 
solubility  of  FeCl2.  In  addition  to  the  solubility  of  iron  chloride,  this 
model  also  predicts  the  movement  of  iron  inside  the  cell  with  time. 
It  has  seen  that  during  discharge  the  solubility  of  FeCl2  decreases 
near  rL  as  xA  decreases.  The  flux  of  the  soluble  ferrous  complex 
becomes  significant  at  moderate  DODs  leading  to  a  non-zero  elec¬ 
trolyte  concentration  of  the  complex  although  solid  iron  chloride  is 
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